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We present a hydrodynamic assessment of preliminary particle spectra observed in Au+Au collisions 
at -y/SNN = 200 GeV. The hadronic part of the underlying equation of state is based on explicit con- 
servation of (measured) particle ratios throughout the resonance gas stage after chemical freezeout 
by employing chemical potentials for stable mesons, nucleons and antinucleons. We find that under 
these conditions the data (in particular the proton spectra) favor a low freezeout temperature of 
around ~ 100 MeV. Furthermore we show that through inclusion of a moderate pre-hydrodynamic 
transverse flow field the shape of the spectra improves with respect to the data. The effect of the 
initial transverse boost on elliptic flow and the freezeout geometry of the system is also discussed. 



I. INTRODUCTION 

During its second year of operation, RHIC (the Rela- 
tivistic Heavy Ion Collider at Brookhaven National Labo- 
ratory) has collided ^^^Au nuclei at center-of-mass {CM) 
energies of 200 GeV per nucleon pair to create strong 
interaction matter at high energy densities in the labo- 
ratory. To identify signals of a possible phase transition 
from low-energy nuclear to deconfined quark-gluon mat- 
ter, a large amount of data was analyzed and recently 
presented for the first time [1] . 

In the present article we investigate single particle 
spectra of various hadronic species within a hydrody- 
namic framework for the reaction dynamics, which as- 
sumes rapid thermalization in the reaction volume and 
a subsequent expansion according to the conservation of 
energy, momentum, entropy and baryon-number (for de- 
tails of the approach, which resides on explicit longitu- 
dinal boost invariance, cf. Ref. [2]). At CM energies of 
130 AGeV, the successful description of observed single- 
particle transverse momentum (px) spectra and their az- 
imuthal modulation in non-central collisions have vali- 
dated this approach down to decoupling temperatures of 
~ 130 MeV, at which hadronic interactions have been 
assumed to cease instantaneously [3]. As an alterna- 
tive to entirely hydrodynamic simulations, especially for 
the late, more dilute stages in a heavy-ion collision, hy- 
brid models have been developed [4-6] which treat the 
hadronic phase in sequential-scattering models, propa- 
gating hadrons individually. While the momentum-space 
observables are in good agreement with experiments at 
RHIC in both descriptions, the freezeout geometry per- 
sists to be inconsistent with the data in either approach. 
Hydrodynamic evolutions appear to be too long-lived but 
too small in radial extent [7], whereas hybrid calcula- 
tions produce an emission cloud which appears to be too 
large [4,5]. 

Concerning global particle production, it was soon re- 
alized [8] that, also at RHIC energies, measured hadron 
ratios reflect a chemical composition of the fireball which 
corresponds to a temperature close to the expected QCD 



phase boundary, T^hem ^ 170-180 MeV ~ T^. Thus, 
in a thermodynamic description of the cooling process 
from chemical to thermal freezeout, the conservation 
of the relative hadronic abundances requires the in- 
troduction of (effective) chemical potentials [9-14] for 
species that are stable on the scale of typical fireball 
lifetimes. In particular, it was pointed out in Ref. [13] 
that the conservation of antibaryons plays an important 
role at collider energies. Despite their large annihila- 
tion cross sections, their finally observed abundance is 
in complete agreement with chemical-freezeout system- 
atics (for a possible microscopic explanation of this fact, 
based on multi-meson fusion reactions to maintain de- 
tailed balance, cf. Ref. [15]). This implies the build- 
up of large antibaryon chemical potentials, , defined 
via /ij^ = ^/^N + M^- Towards thermal freezeout this, 
in turn, entails rather large baryon chemical potentials 
{fiN — 350 MeV), and is at the origin of appreciable pion 
chemical potentials {fi-^ ~ 80-100 MeV). The influence of 
chemical potentials on the hydrodynamic evolution and 
resulting observables has been investigated in Ref. [16] 
for CM energies of 130 AGeV. In addition to conserving 
TT-, K-, rj- and Ty'-numbers, we here explicitly distinguish 
chemical potentials of baryons and antibaryons along the 
lines of Ref. [13] to correctly account for the finite net- 
baryon density at full RHIC energy (200 AGeV). For 
consistency with previous analyses [2,3,17,18] we assume 
a phase transition from quark-gluon to hadron matter at 
Tc = 165 MeV with a latent heat of eiat = 1.15 GeV/fm^ 
and a hadronic resonance gas equation-of-state (EoS) as 
before. At Tc the hadronic phase starts in chemical equi- 
librium to (approximately) reproduce the measured par- 
ticle ratios [8], see above. To improve on previous analy- 
ses, the subsequent hadronic evolution is now constructed 
incorporating effective meson and (anti-) baryon chemi- 
cal potentials as in Ref. [13] to preserve the correct (ab- 
solute) particle abundances. 

As a second new aspect of the present manuscript, we 
present an attempt to refine the initial conditions of the 
hydrodynamic evolution. More specifically, we will ex- 
plore ramifications of pre-equilibrium collective behav- 
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ior by introducing appropriate radial velocity profiles at 
the time of complete thermalization. Such effects can 
be associated with pre-thermal re- interactions, a free- 
streaming period, or a combination thereof, and turn out 
to generally improve the description of transverse mo- 
mentum spectra of the produced particles. 

Our article is organized as follows. In Sects. II and III 
we analyze the impact (and interplay) of off-equilibrium 
hadro-chemistry and modified initial collisions on trans- 
verse momentum spectra of pions, kaons and (anti-) pro- 
tons, both for central and more peripheral collisions in 
comparison to preliminary data at 200 AGeV. Perti- 
nent predictions for azimuthal anisotropics in non-central 
collisions are presented in Sect. IV. We furthermore 
comment on implications for the freezeout geometry in 
Sect. V, and summarize in Sect. VI. 



II. PARTICLE SPECTRA 
COLLISIONS 



CENTRAL 



Let us start by briefly discussing the initial condi- 
tions of our hydrodynamic calculations. According to 
the ~ 15% larger hadron multiplicity at midrapidity in 
central collisions at 200 AGeV [19,20] as compared to 
130 AGeV, we increase the maximum entropy-density 
parameter from Sq = 95 fm^'^ [7] to 110 fm^'^ (keeping 
the equilibration time fixed at tq = 0.6 fm/c to facili- 
tate the interpretation of observed changes). The correct 
baryon admixture is obtained by adjusting the entropy- 
per-baryon to S/B = sq/uq = 250, constant through- 
out the evolution (sq and no are the initial entropy- and 
baryon-density in the center of the collision, S and B the 
total entropy and net baryon number). The thermody- 
namic fields in the transverse plane are set to scale with 
a combination of wounded nucleon and binary collision 
profiles as elaborated in Refs. [7,18], which allows for a 
geometrical prescription to reproduce the multiplicity in 
collisions at finite impact parameter b. 

The results of our calculations with improved hadro- 
chemistry are compared to (preliminary) data for tt^, 
and antiproton px-spectra from central Au+Au col- 
lisions at 200 AGeV [21,22] in Fig. 1 (the experimental 
centrality selection of 5 % is approximated by using an 
average impact parameter b = 2.4 fm). Compared to 
particle spectra in standard (i.e., chemical-equilibrium) 
hydrodynamics we find a better description of the over- 
all curved shape of the hadronic spectra, in particular for 
low-pT pions. This is a result of the meson chemical po- 
tentials (^TT w 80-100 MeV at freezeout), which ampHfy 
the Bose-statistics effect. In addition, the population of 
heavy resonances also increases after inclusion of chem- 
ical potentials which entails larger contributions at low 
Pt from their decay products. At large transverse mo- 
menta the hydrodynamic calculations deviate from the 
data which is suggestive for the onset of the hard scat- 
tering regime. At exactly which values of pt this occurs, 



and how this transition depends on the particle species, 
are among the major questions to be clarified. E.g., high 
energy partons evolving within a hydrodynamic back- 
ground can be introduced to study the particle spectra 
beyond the collective behavior [23]. 

As was already observed in Ref. [16], the expansion 
of the chemically non-equilibrated hadron gas leads to 
slopes for pion spectra that are almost insensitive to 
the decoupling temperature. Proton spectra, on the 
contrary, clearly favor a freezeout at T ~ 100 MeV 
(thick solid line), which corresponds to an energy density 
e ~ 0.075 GeV/fm"^ (which is about the same as in previ- 
ous calculations). The thin lines in Fig. 1 correspond to 
decoupling at the phase transition (recall that the mul- 
tiplicity of the individual particle species is independent 
of freezeout due to the chemical potentials) . 

The experimental pion spectra in the 1-2 GeV range 
appear flatter than what follows from the flow generated 
by hydrodynamic expansion with our given initial con- 
figuration (at transverse momenta pt > 2 GeV this is 
conceivably due to additional perturbative hard scatter- 
ing contributions). To a lesser extent, this is also true for 
the heavier kaons and protons, even at the low freezeout 
temperature of 100 MeV. 
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1. 7r~, K~ and antiproton spectra for central colli- 
sions at 200 AGeV {K~ and p spectra are scaled by factors 
of 1/10 and 1/100, respectively). The thick lines represent 
the results for Tdec = 100 MeV, the thin lines for 165 MeV. 
All calculations are for a thermalization time tq = 0.6 fm/c, 
either without (solid lines) or with (dashed lines) an initial 
transverse boost (see text). 

The data thus seem to exhibit somewhat stronger col- 
lective expansion than developed subsequent to an equi- 
libration time of To = 0.6 fm/c. Additional radial flow 
could be generated by assuming still shorter equilibra- 
tion times, e.g., tq = 0.2 fm/c [24]. It is, however, hard to 
imagine that particles are 'born' into thermal equilibrium 
without allowing for some relaxation time with rescat- 
tering. But even the other extreme, i.e., a period of free 
streaming, induces a non- vanishing radial velocity profile 
due to a separation of originally random particle veloc- 
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ities [25,7]. A realistic situation is probably in between 
the two extremes, essentially pre-equilibrium in character 
with associated rather complicated structures of the gen- 
erated flow-field and energy-density distributions (more 
exotic phenomena such as sphaleron explosions [26] could 
also play a role). As an exploratory study, we here intro- 
duce a simplistic initial 'seed' transverse velocity accord- 
ing to VT{r) — tanh(ar), where r is the radial distance 
from the origin, superimposed on the original fields at 
To = 0.6 fm/c. For a value of a = 0.02 fm~^ the initial 
velocity field for r± < 6 fm/c is similar in magnitude 
(although less parabolic) to both (i) starting the hydro- 
dynamic evolution at earlier time (tq = 0.2 fm/c as in 
Ref. [24]) and evolving it to tq ~ 0.6 fm/c, as well as (m) 
free streaming from r = 0.2 to 0.6 fm/c. The essential 
difference between (i) and (m) lies in the azimuthal dis- 
tribution at Tq = 0.6 fm/c, to which we will come back to 
in Sect. IV. It should also be noted that stronger trans- 
verse flow due to larger transverse pressure is expected if 
the longitudinal expansion is not fully thermalized [27]. 

The results with our simple ansatz are represented by 
the dashed lines in Fig. 1, and are found to improve the 
agreement with experiment, up to pt — 2(3.5) GeV for 
pions and kaons (antiprotons). We note that when in- 
creasing a to 0.05, the proton spectra become much flat- 
ter than experimentally observed. 

III. PARTICLE SPECTRA - NON-CENTRAL 
COLLISIONS 
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FIG. 2. Centrality dependence of positive pion spectra at 
mid-rapidity {y — 0) in terms of central, semi-central (scaled 
by 1/2) and peripheral collisions (scaled by 1/3). 

In Fig. 2 we compare preliminary spectra of positive pi- 
ons [21] to our hydrodynamic results at Tdcc = 100 MeV 
in 3 different centrality bins ('central', b = 2.4 fm, 
iVpart = 343.8; 'semi-central', 6 = 7 fm, A^^part = 170.8; 
'peripheral', b — 9.6 fm, A'part = 76.6). Again, we display 
calculations with an initial transverse boost by dashed 
lines. As expected, the prerequisites for a hydrodynamic 



approach (strong rescattering and a sufficiently large sys- 
tem size) are increasingly invalidated at large impact pa- 
rameters, reflected by an onset of deviations from experi- 
ment at smaller transverse momenta (higher-momentum 
particles can rapidly escape the fireball without thermal- 
izing). For peripheral collisions the agreement between 
theory and experiment holds for pT < 1 GeV, which, nev- 
ertheless, still accounts for more than 96% of the emitted 
particles. 

Fig. 3 shows experimental [21] and calculated proton 
spectra which are of particular interest in the present 
context as they acquire the largest chemical potentials 
{e.g., around thermal freezeout /iN = 380 MeV and 
/iN = 343 MeV implying fi^ = 723 MeV, which yields an 
antiproton-to- proton ratio of 0.72 consistent with experi- 
ment [28]), and are most sensitive to collective expansion. 
We find good agreement of theory and experiment at a 
freezeout temperature of 100 MeV up to pt ~ 3.5 GeV 
in the central, but only up to '--^ 2 GeV in the periph- 
eral sample. The additional transverse 'kick' in the ini- 
tial state as described above (dashed lines) is particularly 
significant for central collisions. 
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FIG. 3. Midrapidity proton spectra for central, semi-cent- 
ral (scaled by 1/2) and peripheral collisions (scaled by 1/3). 

Despite the fact that particle densities in the later 
hadronic stage of the expansion are moderate, we con- 
clude that rescattering is strong enough to allow for a 
hydrodynamic description until thermal decoupling. The 
correct chemical composition of the hadronic gas is main- 
tained by the generation of large chemical potentials, 
which (at given temperature) provide an increased num- 
ber of scattering partners with larger cross sections as 
compared to a chemically equilibrated environment. 

We have focused here on positively charged pions and 
protons. The corresponding results for tt", K^, and 
p in non-central collisions are of similar quality. The mul- 
tiplicities and mean transverse momenta of these parti- 
cles are collected in Table 1. 
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rv — n on 
1 2 3 


— U.U^ ill! 

1 2 3 


dN 

(pt) 


280.8 134.3 57.84 
0.398 0.392 0.375 


282.3 134.6 57.82 
0.419 0.405 0.383 


dN 

K+ dy 

(pt) 


50.18 23.99 10.33 
0.619 0.608 0.572 


50.43 24.05 10.33 
0.660 0.634 0.589 


dN 
p dy 

(pt) 


28.08 13.44 5.798 
0.880 0.861 0.802 


28.13 13.44 5.794 
0.949 0.906 0.831 



Table 1; Multiplicities and mean transverse momenta (in 
GeV) of different particles for 3 centrality selections (l=cen- 
tral, 3=periplieral) at y = 0. The p/p ratio is 0.72; (pt) of 
antiprotons is within 1% of the proton value. 



IV. ELLIPTIC FLOW 

For the same impact parameters as considered above 
we proceed by studying the azimuthal anisotropies of par- 
ticle spectra [29], i.e., the momentum dependence of el- 
Uptic flow as defined by W2(pt; b) = (cos(20)), where the 
average is taken over the angular distribution of particles, 
dN/ dypTdpTdcj). 

Flow anisotropy is generated during the earliest stages 
of the collision, at which the spatial eccentricity of the 
thermodynamic fields and the anisotropies in the pres- 
sure gradients are the largest. The matter is set into 
anisotropic motion as larger forces are acting along the 
'short' radius of the initial (overlap) ellipse. This motion 
rapidly reduces the spatial anisotropies, thereby bring- 
ing further generation of momentum anisotropy (i.e., 
V2) to a stall [30,2]. If the system evolves in chemi- 
cal equilibrium, the dominant particle species at freeze- 
out are pions, which carry the generated anisotropy in 
their momentum distribution. Their differential ellip- 
tic flow, W2(pt), is then almost independent of the de- 
coupling temperature Tdoc [3]. Heavier particles, on the 
other hand, do exhibit some dependence on Tdec, mainly 
because of the continuously increasing radial flow which 
shifts the generated anisotropy towards larger transverse 
momenta. In the presence of effective chemical poten- 
tials the contribution of protons to the total anisotropic 
flow (of all particles) is still small; however, the contri- 
bution of their number to the particle yield is more sig- 
nificant. The anisotropic flow must thus be absorbed by 
the pions (which, due to their small masses, adjust their 
momentum distribution easier) . Through this effect their 
elliptic flow now also becomes sensitive to the decoupling 
temperature, as found in Ref. [16]. In addition, the influ- 
ence of resonance decays is enhanced in the chemical off- 
equilibrium formulation. The heavy resonances, which 
at large transverse momentum carry rather large ellip- 
tic flow, decay and transfer their elliptic flow to pions at 
relatively low transverse momentum. 



In Fig. 4 we show results for elliptic flow of pions 
(left panel) and protons (right panel) from the hydrody- 
namic calculation under inclusion of chemical potentials. 
The initial transverse boost as defined in Sect. II shifts 
the anisotropy to larger transverse momenta which im- 
plies a reduction of V2 at given p^. The development of 
anisotropic fiow is additionally hindered since it has to 
form on top of the isotropic initial boost field which we 
have employed here. The value Cp at which the anisotropy 
of the energy-momentum tensor of the fluid T^^" satu- 
rates during the evolution [2] , is about 25% smaller than 
without the initial 'kick'. 

It is instructive to compare these results to the sce- 
nario where the equilibration time is set to very small 
values (tq = 0.2 fm/c for the dashed-dotted curves in 
Fig. 4). As elucidated in Sect. II, this generates as much 
radial flow as the superimposed profile at tq = 0.6 fm/c 
does. However, the elliptic fiow is larger than in the for- 
mer case, but not significantly different from using an 
equilibration time of 0.6 fm/c without initial kick. This 
is due to the fact that, without initial kick, V2 saturates 
for either equilibration time at approximately the same 
value *. 
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FIG. 4. Elliptic flow of positively charged pions (left) and 
protons (right) for three different impact parameters. The 
dashed lines include an initial transverse boost as described 
in the text. Dashed-dotted lines represent the results when 
assuming thermalization at ro = 0.2 fm/c with a = 0. 

The experimentally observed elliptic fiow reaches a lim- 
iting maximal value as a function of transverse momen- 
tum. The PHENIX collaboration has pointed out [21] 
that this saturation is reached at smaller transverse mo- 
menta for pions than for protons, and that the satura- 
tion value appears to be larger for the latter. Within the 
hydrodynamic framework this refiects the earlier break- 
down of the strong rescattering assumption for pions, 
which for protons remains valid up to higher pT due to 
larger (average) scattering cross sections (o-TrAf > (i^Tr). 
This is corroborated by the description of the single- 



*Note that this is no longer true for significantly larger tq, 
e.g. 2 fm/c, for which V2 is significantly reduced and under- 
predicts the data already at 130 AGeV 
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particle spectra which extends to larger for (anti-) pro- 
tons than for pions. The different z;2-saturation momenta 
for mcsonic and baryonic elliptic flow are also consistent 
with the formation of hadrons via quark-coalescence [31]. 
Within this picture one similarly expects a larger satu- 
ration value of V2 for protons than for pions. 

V. FREEZE-OUT GEOMETRY 

Let us finally comment on the implications of our re- 
sults for the frcezeout geometry of the hadronic system. 
In Ref. [14] it was pointed out that the relation between 
energy density and pressure, e{p), for the hadronic equa- 
tion of state is barely modified by the introduction of 
chemical potentials. Therefore, the space-time evolution 
of the system, which is largely driven by this relation, is 
not substantially altered either. A large change, however, 
occurs in the relation between temperature and energy 
density, T'(e), which thus influences the construction of 
the freezeout hypersurface and the thermal properties of 
the fluid on this surface. E.g., in chemical equilibrium 
the energy density at T = 130 MeV corresponds to a 
temperature of only 100 MeV in the presence of large 
chemical potentials, since the latter increase particle and 
energy densities approximately by pertinent fugacity fac- 
tors e^/'^. Therefore, the freezeout hypersurface of the 
hydrodynamic calculations in chemical off-equilibrium is 
not much different from the hypersurface of previous cal- 
culations if freezeout is performed at a comparable energy 
density (i.e., the freezeout temperature is adapted ac- 
cordingly). In both cases, the fireball decouples at about 
15 fm/c after equilibration (in central collisions) and has 
about the same spatial extent. Only after inclusion of 
the initial radial flow proflle is the lifetime shortened by 
~ 15%, and the transverse expansion increases by about 
the same percentage. For obscrvables, this entails smaller 
longitudinal correlation radii (which reflect the system's 
lifetime) but only slightly larger sideward radii. This 
effect reduces the discrepancies between calculated and 
measured Hanbury-Brown and Twiss (HBT) radii by a 
few percent [7], but is not sufficient by itself. Additional 
eff'ects, such as viscosity [32], large partonic cross sec- 
tions in the early phases [33], or a refined treatment of 
hadronic rcscattering [5] and freezeout [34] (including, 
e.g., a large p-meson width as predicted in Ref. [35]), 
seem to be required to fully resolve the "HBT puzzle". 

VI. SUMMARY 

Based on a resonance gas equation of state which ex- 
plicitly incorporates hadrochemical freezeout by employ- 
ing chemical potentials for (stable) mesons and baryons 
in the hadronic evolution, we have performed hydrody- 
namic simulations of heavy ion-collisions at full RHIC en- 
ergy. We have compared the results for pion, kaon, and 



proton j?T-spcctra to preliminary data from 200 AGeV 
Au-|-Au collisions at different centralities. Our in- 
vestigations indicate the necessity of an initial (pre- 
hydrodynamic) transverse flow to better account for the 
slopes of the observed spectra. Good agreement with pre- 
liminary data for transverse momentum spectra in cen- 
tral collisions is obtained up to ~ 1.5-2 GeV for pions, 
and up to at least 3 GeV for protons. We further studied 
the influence of hadrochcmistry and initial flow on elliptic 
flow and source geometry. The former has been presented 
as a prediction for pions and protons for upcoming ex- 
perimental analyses. For the latter, some improvement 
with respect to the discrepancy between model and data 
has been found, but additional effects remain mandatory. 
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